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Abstract 



^—^ ■Numerical simulations of dispersive turbulence in magnetized plasmas based on the Hall-MHD description are presented, 
assuming spatial variations along a unique direction making a prescribed angle with the ambient magnetic field. Main 
observations concern the energy transfers among the different scales and the various types of MHD waves, together with 
'the conditions for the establishment of pressure-balanced structures. For parallel propagation, Alfven-wave transfer to 
Q^small scales is strongly inhibited and rather feeds magnetosonic modes, unless the effect of dispersion is strong enough 
at the energy injection scale. In oblique directions, the dominantly compressible character of the turbulence is pointed 
out with, for quasi-transverse propagation, the presence of conspicuous kinetic Alfven waves. Preliminary simulations 
of a Landau fluid model incorporating relevant linear kinetic effects reveal the development of a significant plasma 
temperature anisotropy leading to recurrent instabilities. 
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1. Introduction 

Turbulence in magnetized plasmas remains a main is- 
sue in the understanding of the dynamics of media such 
as the solar corona, the interstellar medium, the so- 
lar wind or the planet magnetosheaths. In the so- 
iar wind for example the turbulent cascade extends 
much beyond the ion Larmor radius. One of the ques- 
tions concerns the spectrum of the magnetic fluctua- 
tions that displays a power-law behavior on a broad 
jange of wave numbers, with a co nspicuous change of 



slope near the inverse io n gyro r adius ([Leamon et c 



bolstein and Roberts!. 1999t Alexandrova et al 



1998 



2006 



Sahraoui et al.l . 120091 ) . This effect is often associated with 



.the i nfluence of wa.ve dispersion , induced by the Hall cur- 
rent ( Ghosh et al. ■ 11996 1: iGaltieil I 20061: Alexandrova et al 



Shaikh and Shuklal . 
superposition of 



!2007t ICaltier and Buch hn. 2007 t IServidio et al.L l2'o07 : 
^009|), but could also result from 
cascades of kinetic Alfven waves 
and ion entropy fluctuations, as suggested by studi es 



based on the gyrokinetic f ormalism (|Howes et al 



Schekochihin et al.l . [2009h 



At scales large compared with the ion inertial length or 
the ion Larmor radius, the usual MHD description pro- 
vides a satisfactory description of regimes where, due to 
the presence of a strong ambient field, a dominant effect is 
the anisotropic energy transfer to Fourier modes with large 



Email addresses: lavederSoca. eu (D. Laveder), 
marradiSoca.eu (L. Marradi), passotOoca. eu (T. Passot), 
sulemSoca.eu (P.L. Sulem) 



transv er se wavenumbe rs (see e.g. iGhosh and Golstein 



( 19971) : lOughto n and Matthae"ua (|2005l ) and references 
therein). This suggests that the dynamics of transverse 
small sc ales may b e amenab le to a reduced MHD descrip- 
tion ('(Z ank and Matthaeus . 19921) and references ther ein), 
possibly including Hall current ( Gomez et al.l . l2008l) or, 
when retaining scales significantly smaller th an the ion 
Larmor radius, to a gyrokinet ic approach ((Howes et al 



2006t ISchekochihin et all . l2009l ). The latter that appears 



to be very efficient in describing strongly-magnetized near- 
equilibrium fusion plasmas is still under discussion con- 
cerning its applicability t o space and astrophysical plasmas 



( Matthaeus et al. . 20081) . In the solar wind for example 



magnetic fluctuations may be comparable to the ambient 
field. Furthermore, longitudinal transfer could a priori be 
non negligible in a compressible regime, at scales where 
Hall current and kinetic effects play a significant role. A 
weak turbulence theory performe d on the Vlasov-Maxwe ll 
system was recently developed ( Yoon and Fanj . 12008 ). 
showing the existence of a parallel cascade of low-frequency 
Alfven waves through a three- wave decay process mediated 
by ion-sound turbulence, in a regime where wave-particle 
interactions are neglected. Addressing this issue by di- 
rect numerical simulations of the Vlasov-Maxwell equa- 
tions being still difficult on the present-day computers, the 
question arises whether a similar cascade can be observed 
within a fiuid model that retains important ingredients 
of the above theory, such as compressibility and disper- 
sion. As a first step, we address the problem within the 
simplest description provided by Hall-MHD (HMHD) with 
Ohmic and viscous dissipations, together with a large-scale 
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external driving acting on the transverse components of 
the velocity or magnetic field. We specifically concentrate 
on a one-dimensional setting where the variations of the 
fields are restricted to a direction making a prescribed an- 
gle with the ambient magnetic field, a framework that al- 
ready reveals a manifold of complex dynamical processes 
that deserve detailed investigations before including addi- 
tional physical and multidimensional effects. In the case 
of quasi-transverse propagation, we also present simula- 
tions of a model that extends the HMHD by retaining 
pressure anisotropy, Landau damping and finite Larmor 
radius effects up to transverse scales significantly smaller 
th an the ion Larmo r radi us. This approach developed 
in iPassot and SulemI (2007) extends the so-called Landau 



fluid model initiated in ISnvder et all (|f997f ) for the MHD 
scales where Landau damping is the only relevant kinetic 
effect. 

The paper is organized as follows. Section 2 briefly re- 
views the Hall-MHD description and its one-dimensional 
reduction. Section 3 concentrates on the case where the 
dynamics takes place in the direction of the ambient field. 
The case of oblique propagation is addressed in Section 
4. Landau fluid simulations retaining small-scale kinetic 
effects are reported in Section 5. Our conclusions are sum- 
marized in Section 6. 

2. The Hall-MHD description 

HMHD can be viewed as a bi-fluid description of a 
plasma, where electron inertia is neglected. The presence 
of the Hall term in the generalized Ohm's law allows a 
decoupling of the ion fluid from the electron one in which 
the magnetic field lines are frozen. The va lidity conditions 
of HMHD arc discussed in Howed (j2009(l where compar- 
isons with kinetic theory are presented. Choosing as units 
the Alfven speed, the amplitude of the ambient magnetic 
field, the equilibrium density and the ion inertial length 
li (defined as the ratio of the Alfven speed to the ion gy- 
rofrequency) , the HMHD equations (for the ion fluid) read 



9tp + V.(pv) = 



(1) 



f3 



p(9tv + v.Vv) = -^V/9'' + (V X b) X b (2) 
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^tb - V X (v X b) 
V.b = 0, 



-V X (^-(V X b) X bj(3) 
(4) 



where the total /? parameter is the square ratio of the 
sound to Alfven velocities, and a polytropic equation of 
state p (X is assumed for both ions and electrons. 

When the spatial variation is restricted to a dependency 
on the X coordinate along a direction making an angle 9 
with the ambient magnetic field Bo = (cos 6*, sin 6*, 0), one 
gets 

dtp + dx(pv^) ^ (5) 
dtv, + v^d^v, = --dJ^p^ + + bl)) 



p 



P 



(6) 



(7) 
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± cos Odrr I -dxb[z^y] \ + K[y^z-\dxxb[y,z\ 



fl. 



where driving and dissipation have been supplemented in 
both the velocity and magnetic field equations. Here, the 
subscript [y, z] refers to the vector component along the y 
or the z direction, or to the value of the viscosity acting 
on the corresponding velocity component. The ± sign in 
front of the Hall term depends on the considered compo- 
nent of the magnetic field. No (artificial) hyperviscosity 
and magnetic diffusivity nor spectral filtering are used in 
the simulations. Instead, anisotropic dissipations are as- 
sumed. In the case of parallel propagation, different vis- 
cosities and diffusivities are taken in the directions par- 
allel and transverse to the ambient field, by prescribing 
Ky = Kz — fJ-y — P'z fJ-x- For oblique propagation, we as- 
sume smaller coefficients in the direction perpendicular to 
the plane defined by the magnetic field and the direction 
of propagation, in the form Hz — Pz Px = Py = i^y 

The driving is assumed to act either on the velocity 
(kinetic driving) or the magnetic field (magnetic driv- 
ing) components. For parallel propagation, we prescribe 
fy = fz oi' fy = fz while, for oblique propagation, the 
driving reduces to or to Such a driving is supposed 
to minimize the sonic components, as it is acting on field 
components perpendicular to the ambient field in parallel 
propagation and to the plane defined by the ambient field 
and the propagation direction when the latter is oblique. 
The values of the diffusivity and viscosities in the various 
directions are chosen as the minimal values (depending on 
the spatial resolution and of the physical parameters of the 
runs) needed to accurately resolve all the retained scales. 

In all the simulations, we take 7 = 5/3 and (3 — 2. Each 
component of the kinetic or magnetic driving (generically 
denoted /) is a white noise in time defined by its Fourier 
transform fk — CS, ^Fk/At where ^ is a Gaussian random 
variable with zero mean and unit variance, chosen indepen- 
dently at each time step. This ensures a constant mean 
flux of energy injection that can be chosen at will, as in the 
usual phenomenology of the turbulent cascades. Further- 
more, such a driving process avoids an artificial enhance- 
ment of a specific type of waves and enables the emergence 
of the dominant modes as the result of the nonlinear dy- 
namics. The spectral distribution i^fe = k* exp(—{2k'^ / k'j)) 
is peaked about a wavenumber kf. 

The HMHD system is integrated in a periodic domain 
using a Fourier pseudo-spectral method where most of the 
aliasing is removed by spectral truncation of the computed 
nonlinear terms at 2/3 of the maximal wavenumber. The 
spatial resolutions given in the following sections are the 
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effective ones, after aliasing has been suppressed. In all 
the simulations, the te mporal scheme is a. third-order low- 
storage Runge-Kutta ( Williamson . 1980l) . Resolving all 
the temporal scales present in the system, this scheme ac- 
curately preserves the dispersion relation of all the linear 
modes retained in the sim ulation, in contrast w ith implicit 
or semi-implicit schemes ( Laveder et al. . 20091 ). 



For convenience, we collected in Table 1, the main pa- 
rameters characterizing the simulations discussed in the 
forthcoming sections. 

As seen in the following, in spite of the turbulent regime 
achieved in the HMHD simulations discussed in this paper, 
signatures of the linear waves are often present. It is thus 
useful to briefly review the linear theory of eigenmodes 
for the HMHD equations in the absence of dissipation and 
driving. By linearizing eqs. ([S])-® about the equilibrium 
state associated to p = 1, bx = cos 9, by — sin 9, one derives 
that the (real) eigenfrequencies uji with i — 1,2,3 obey the 
dispersion relation 

- fc2(l + /3 + cos^ 9 + k^ cos^ 9)uj^ 



the LH polarized waves become ion-cyclotron waves. It 
is noticeable that the dispersive parallel Alfven waves are 
also exact solutions of the nonlinear HMHD problem. 

For an oblique direction of propagation and /3 > 1, we 
still refer to the eigenmode corresponding to the smallest 
eigenvalue oji, that is dominated by the field components 
Vy, Vz (fig. [H, as the ion-cyclotron wave. In the case of 
quasi-transverse propagation, it is often called the kinetic 
Alfven wave. The intermediate mode corresponding to 
the eigenvalue uj2 is the one that reduces to a pure sonic 
mode for 6* = 0. For oblique propagation it is dominated 
by the component at small scales only, while at large 
scales bz and Vz prevail. The mode at highest frequency 
wa, which corresponds to the whistler wave for 9 = 0, has 
properties opposite to those of the intermediate one. For 
oblique propagation, it is usually viewed as a whi stler wave 



at small scales only (jKrauss-Varban et al.l . 119941 ) . For the 



sake of simplicity, we will keep this terminology at large 
scales also, where this mode is strongly compressible. 



cos^ 61(2/3 + 1 + /3fc2)tj2 _ ^fc6 ^,Qg4 Q^Q^ 3. Parallel dynamics 



where it is sufficient to concentrate on the positive solu- 
tions. 

Assuming an oblique propagation, the corresponding 
eigenmodes are given by 



W(^) = (p, vx = a^lp, Vy = a«p, Vz = a«p. 



= <^blP^ bz 



{i) \ 
O^b P) 



with 



a), 



1 



(i) 



(fc^ cos^ 9 — Lof) sin 9 
k^cos'^9- (^^cos'^9~Lo 
—ik'^Ldi sin 9 cos 9 



fc4 cos2 9 



(S 



cos^ 9 — LJi 



cos 9- 



-a 



a 



(i) 



coa9- 



k 



-a 



(10) 

(11) 

(12) 

(13) 

(14) 
(15) 



The variation of these coefficients with the wavenumber k 
is displayed in fig. [Ufor 9 = 45° and 9 = 80°. 

For parallel propagation {9 = 0), the sonic wave for 
which uj = \/]3k decouples from the eigenmodes corre- 
sponding to (circularly-polarized) Alfven waves. The lat- 
ter obey by + ibz = —{uj/k){vy + ivz) cc eiq)[ia{kx — ujt)] 
where cr = —1 or 1, depending of the left-hand (LH) or 
right-hand (RH) polarization, with the dispersion relation 
= (tT/2)fc^ + fc(l + fc^/4)^/^. In the large wavenumber 
limit, the frequency of LH polarized waves saturates, while 
it grows like fc^ in the case of RH polarization. At fre- 
quencies beyond the HMHD asymptotics, the RH polar- 
ized Alfven waves are continued into whistler waves, while 



3.1. Large-scale driving 

We first consider a domain of extension L = IGtt in a 
direction parallel to the ambient magnetic field. We use 
a transverse kinetic or a magnetic driving characterized 
by the parameters C = 0.1 and kf = 0.5 correspond- 
ing to a mode index rif = A. Defining the ion inertial 
wavenumber as ki = 27: /U, this corresponds to a ratio 
ki/kf = An. A spatial resolution of 1024 grid points is 
used, which prescribes a time step as small as At = 5-10"^, 
due to the dispersion relation of the whistler modes. We 
use a viscosity p.^ = 10^^ in the equation for the paral- 
lel velocity, and equal viscosity and magnetic diffusivity 
IJ'y = iJ'z = = K,y = 10~^ for the trausvcrsc fields. 
This relatively large viscosity in the parallel direction is 
required because of the development of shock waves. 

In the case of kinetic driving (run A), fig. [2] (left) 
shows the time evolution of the total energy of the sys- 
tem (from which the initial value has been subtracted) 
E = /(! |vp + i(|bp - 1) + ;^(p^ - l))dx. It 
also displays the contributions of the transverse kinetic 
EX_ = /f|v_Lpd2: and magnetic E'^ = ^ \\\n_\'^dx en- 
ergies that turn out to be comparable and much larger 
than the parallel kinetic and the internal energies (not 
shown). It is conspicuous that the energy does not sat- 
urate, even after an integration time t = 20000. Figure 
[2] (right) displays the time evolution of the magnetic en- 
ergy on the modes of index 1, 2, 3, 4 (a similar behavior is 
obtained for the kinetic modes). We observe that after a 
while, the mode of index n = 1 becomes strongly domi- 
nant, indicating a significant inverse transfer to larger and 
larger scales (leading to a non saturation of the energy), as 
expected from an equ ilibrium thermodynamical argument 
( Servidio et al.l . I2OO8 ) , up to the moment when it reaches 
the size of the computational domain. 
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run 


domain size 


propagation angle 


driving 


A 


L = IBtt 


61 = 0° 


kinetic, C = 0.1, kjU = 1/2 


B 


L = 167r 


= 0° 


magnetic, G = 0.1, = 1/2 


C 


L = 47r 


61 = 0° 


kinetic, C = 6.25 x 10"^, kfk = 2 


D 


L = 47r 


61 = 0° 


magnetic, C = 6.25 x 10"'^, k-fk — 2 


E 


L = 167r 


e = 45° 


kinetic, C = 0.1, kfk = 1/2 


F 


L = 47r 


e = 45° 


kinetic, C = 6.25 x lO'^, kfk = 2 


G 


L = Att 


61 = 45° 


magnetic, C = 6.25 x 10"^, kfk = 2 


H 


L = l&TT 


e = 80° 


kinetic, C = 0.1, kfk = 1/2 



Table 1: Simulation parameters of HMHD simulations 



A main observation is the establishment of a pressure- 
balanced state (fig. [3j left) that is not limited to the largest 
scales, in the sense that pressure balance is still observed 
when the largest Fourier modes are filtered out. Inspec- 
tion of the square transverse velocity |v_Lp (fig. [31 right) 
shows the presence of small scales superimposed to the 
n — 1 mode, which are more conspicuous than on the mag- 
netic component \h±\'^ displayed on the left panel. When 
filtering out the modes of index n = 1, the square trans- 
verse velocity |vj_p reveals the presence of solitonic waves 
which in many instances survive collisions and preserve 
their coherence on several time units (fig. IH left). It is of 
interest to notice that dispersive pressure-balanced struc- 
tures are commonly ob served in space plasmas (see e.g. 

20031 )). The right panel that displays 



(Stasiewicz et al 



the individual transverse velocity components also shows 
the presence of rotational discontinuities (near x = Air), 
together with the existence of fluctuations at very small 
scales. 

Figure [S] displays the transverse kinetic |vj,p + jv^p 
(thick line) and magnetic jb^p -f jb^p (thin line) spectra, 
averaged over the time intervals t — 18200 — 18250 (left), 
and t — 19500— 19550 (right). The left panel corresponds 
to a stage of the evolution for which flg. H] provides instan- 
taneous snapshots. We can distinguish various spectral 
ranges, and in particular two different power-law domains 
at large and intermediate scales whose exponents fluctu- 
ate in time while preserving a kinetic spectrum shallower 
than the magnetic one. Such a dominance of the kinetic on 
the magnetic contribution suggests an ion-cyclotron turbu- 
lence since for these waves the transverse velocity and mag- 
netic field components are in a ratio k/uj > 1. The transi- 
tion between the two spectral ranges occurs near k ^ ki = 
2tt (reminisc ent of observati o ns in t he solar wind and mag- 
netos h eath ([Leamon et aL , 1998 ; Golstein and Roberts! . 



1999t lAlexandrova et all . i200a ISahraoui et all . 1200^ ;) 
that corresponds to the typical scale of the solitonic struc- 
tures seen in fig. |31 The flat region visible at larger 
wavenumbers is the spectral signature of the small-scale 
fluctuations visible in physical space. At longer times 
(right panel), the intermediate power-law range has ex- 
tended as the soliton amplitude increased, taking over the 
flat spectral region. 



In the case of a magnetic driving (run B), a similar over- 
all dynamics is obtained but on a longer time scales (not 
shown). The energy accumulation on the largest scale 
is slower by a factor four to five, and the pressure bal- 
ance is less conspicuous even at the end of the simulation 
{t = 10000). Furthermore, no clear power law spectra have 
developed. 

Further insight on the dynamics is provided by the dis- 
sipation of the various fields. Figure [5] (left) displays 
the time evolution (averaged over a time interval of 1000 
time units) of the viscous dissipation D\\ = J ^xipxVxYdx 
originating from the parallel velocity for the two types 
of driving. The dissipation observed is slightly larger in 
the case of magnetic driving and reflects a stronger di- 
rect energy transfer to sonic waves. This observation is 
consistent with the better pressure balance obtained with 
kinetic than with magnetic forcing. On the right panel 
is displayed the (similarly averaged) dissipation D±^ = 
J iHyidxVy)'^ + fizidxVz)'^ + Ky{dxby)'^ + K^{dxb:^)'^)dx of the 
transverse flelds, that turns out to be larger in the case of 
kinetic driving, indicating a more efficient transfer to small 
scales for the transverse fields. 

Energy dissipation affects dominantly the parallel ve- 
locity component, the transverse dissipation being about 
two orders of magnitude smaller. This suggests that the 
injected energy is mainly transferred to sonic waves and 
dissipated through a cascade of acoustic waves. For par- 
allel propagation, writing the equations for the energy 
density of the Alfven and acoustic waves in the form 

a(fKP + i|bj|) 



1 



bA_\'^dxVx + M_LV_L • dxxV± + K±b_L • dxx^A 



• fl + bi • f 

^ 

7(7- 



(16) 



7 , P 2 

JT'' + 2"- 



-\b^\'^dxVx + PxVxdxxVx + Vxfx + bxfx (17) 



where 



hx — —£jlqbjbxdxbq 



(18) 



arises from the Hall term, we get that the energy transfer 
from the Alfven to sonic waves is given by 



\h_i\'^dxVxdx. 



(19) 



Figure [7] displays in the case of kinetic driving the time 
variation of S, together with that of the parallel (viscous) 
dissipation. Left panel, which corresponds to instanta- 
neous quantities, shows that the amplitude fluctuations of 
the transfer are much larger than those of the dissipation. 
We nevertheless observe on the right panel which displays 
the same quantities averaged on a time interval of 1000 
time units, that the parallel dissipation identifies in the 
mean with the energy transfer. A similar behavior is visi- 
ble with the magnetic driving. 

The above observations indicate that sonic wave turbu- 
lence is the dominant phenomenon, although small scales 
also form on the transverse field components on a longer 
time scale. In order to address more precisely the nature 
of the transfers, it is of interest to perform a simulation 
where the final time t = 20000 of the previous simulation 
is taken as initial condition and the energy is no longer in- 
jected at a constant rate, but the driving is monitored in 
order to maintain the total energy almost constant. This 
constraint leads to a drastic change on the transverse spec- 
tra which, within a few thousands of time units, loose en- 
ergy at the scale of the solitons displayed in fig. 21 leading 
to a conspicuous spectral gap (fig. [51 left). In physical 
space, solitonic structures have indeed disappeared, the 
profile of |v^p reducing to small-scale oscillations super- 
imposed to large-scale quasi-sinusoidal structures, as seen 
in fig. [5] (right) where the strongly dominant mode n = 1 
has been filtered out. At scales larger than the spectral 
gap, the pressure-balanced state persists, while at smaller 
scales viscous and Ohmic dissipations are negligible and 
dispersion is significant, thus permitting the small scales 
to persist. This suggests that, in the present simulation, 
the transfer of Alfvenic energy from large to intermediate 
Alfven scales which in the previous simulation maintained 
the solitons, does not take place in the absence of a suf- 
ficiently strong driving, and thus cannot compensate the 
transfer of energy from the Alfvenic to the sonic modes 
at the soliton scales. The dynamics thus turns out to be 
significantly different from ordinary turbulence and should 
rather be viewed as a structure-dominated regime display- 
ing significantly less universal properties. 

A further characterization of the dynamics is provided 
by a frequency analysis on the spatial Fourier modes of 
the fields components. In contrast with the incompressible 
regime where the presence of the sole Alfven waves makes 
convenient the analysis of time records a t given points in 
physical space (jDmitruk and Matthaeud . l2009t ). the pos- 
sible dominance of different modes at different scales in 



the present case, leads to discriminate between different 
Fourier modes. Figure [3] displays the temporal spectrum 
of the spatial Fourier modes k = 1, k — 4 and fc = 25 
respectively, for the components Vx (blue) , bz (red) and Vz 
(green) for kinetic (run A: panels a,b,c) and magnetic (run 
B: panels d,e,f) driving. The analysis has been performed 
during a time interval close to the end of the simulations. 
At large scale {k ~ 1 = 2fc/), two peaks are visible. They 
correspond, up to a slight infrared shift possibly due to 
the presence of large-scale coherent structures, to the ion- 
cyclotron {uj « 0.62) and the sonic {u w 1.4) frequencies. 
These peaks are wider with kinetic than with magnetic 
forcing, indicating stronger nonlinear couplings in the for- 
mer case which is also characterized by a larger amount 
of transverse energy at all scales, with a ratio between the 
kinetic and magnetic components in qualitative agreement 
with the linear ion-cyclotron eigenmode. At smaller scales 
(fc = 25 « 4fci), the relative importance of the sonic peak 
decreases in the case of kinetic driving, but remains sig- 
nificant for magnetic driving, indicating a higher level of 
sonic turbulence in the latter regime. 

3.2. Driving at smaller scales 

In a second series of simulations, the global effect of 
the Hall term has been increased by injecting energy at 
smaller scales, with the aim to investigate the influence of 
dispersion on the energy transfers. To address this issue, 
the extension of the computational domain was reduced 
to L = 47r and the system driven at kf = 2, leading to 
a ratio ki/kf = n. The rate of energy injection was low- 
ered by the same factor as the domain extension by taking 
C = 6.25-10~'^, in order to preserve the same amount of en- 
ergy injection per unit length. Reducing the resolution to 
256 mesh points, ensures the same maximal wavenumber 
as in the simulations of Section 3.1. To preserve numeri- 
cal stability, we had to use time steps At = 3.125 • 10~^ 
and At — 2.5 • 10^^ for the kinetic and magnetic drivings 
respectively. The required viscosities and magnetic diffu- 



10- 



10- 



sivities are now fj,x 
both cases. 

Kinetic (run C) and magnetic (run D) drivings lead now 
to significantly different dynamics. In the former case, the 
inverse cascade is less efficient than in the conditions of 
Section 3.1. The energy of the modes with n ~ 1 still in- 
creases up to the end of the simulation (not shown) , but 
after a time of about 7000 times units, the accuracy of the 
simulation deteriorates, the temporal resolution becoming 
insufficient. The total energy grows slowly but in contrast 
with the conditions of Section 3.1, the transverse magnetic 
energy is now significantly lower than the kinetic one (fig. 
[TUl left). Differently, in the case of a magnetic driving, ki- 
netic and magnetic energies are comparable, and the sys- 
tem reaches a stationary steady state where the energies 
saturate (Fig. [TOl right). In the case of a kinetic forcing, 
pressure-balanced structures are still present, while with a 
magnetic driving, no pressure balance establishes. 
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Figure 9: Frequency spectra (red), (green), \vy\'^(k,Lo) (black), \va:\'^ {k,Lo) (blue) as functions oiu), at various wavenum- 

bers k for the runs described in Sections 3 and 4 (see Table 1). The vertical lines refer to the three linear eigenfrequencies given by the 
dispersion relation (|9]l. When only two lines are visible, they correspond to the two higher frequencies, while when only one is present, it is 
the intermediate one. 
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Figure [Til displays the profiles of \v±\'^ for kinetic (left 4. Oblique dynamics 



panel) and magnetic (right panel) drivings, showing that 
more small scales develop in the former case. The struc- 
tures are traveling much more slowly than with the large- 
scale forcing. Furthermore, when comparing transverse 
velocity and magnetic fields, we note that the latter is sig- 
nificantly smoother for both types of drivings. 

The contrast between the two types of forcing is even 
more striking at the level of the dissipations (fig. [T^ . 
While for magnetic driving, dissipation is strongly domi- 
nant on the parallel magnetic field fiuctuations, with the 
kinetic one it mostly affects the transverse components, 
indicating a reduced transfer to sonic waves and the possi- 
bility of a direct cascade of Alfvenic modes, in qualitative 
agr eement with th e weak -turbulence analysis performed 
by lYoon and Fan3 ( 2008 ) on the Vlasov equation. Fur- 
thermore, when the kinetic forcing is monitored in order 
to maintain a prescribed energy as in Section 3.1, no en- 
ergy gap establishes, indicating that the structures can 
sustain a much more constant direct transfer than in the 
case of Section 3.1. This suggests the possibility of a more 
standard turbulence when the injection scale is sufficiently 
close to that at which structures form and dispersion acts 
efficiently. 

The significant difference between the two forcings can 
be clarified by looking at the distribution of energy among 
the various modes and the different scales, presented in 
fig. [HI At large scale (fc — 4, panel m), a peak is visi- 
ble at the ion cyclotron frequency for both drivings, while 
whistler modes are also present in the case of magnetic 
driving (panel o). The relative intensity of the magnetic 
versus kinetic components is consistent with the proper- 
ties of the corresponding eigenmodes. Note that with a 
magnetic driving, the peaks are especially sharp for the 
transverse modes. When the wavenumber is increased 
{k = 25), the peaks broaden under the effect of wave cou- 
pling. This leads to dominant ion cyclotron turbulence in 
the case of a kinetic forcing (panel n). With a magnetic 
driving (panel p), the sonic mode, visible on the left of the 
whistler peak at large scale, becomes dominant, while the 
whistler peak still remains conspicuous. To summarize, 
the pressure balance observed in the case of a kinetic driv- 
ing is consistent with the small amplitude of the interme- 
diate modes, while the transverse dynamics is dominated 
by ion-cyclotron modes that, with the present forcing, can 
cascade down to the dissipation scales. Differently, in the 
case of a magnetic driving, the transverse velocity is dom- 
inated by ion-cyclotron modes and the transverse mag- 
netic field to whistler modes. The large amplitude of the 
latter waves prevents in particular the establishment of a 
pressure balance when the system is magnetically driven, 
which leads to an efficient generation of sonic modes that 
cascade and dissipate. In contrast, kinetic forcing does not 
drive whistler waves, so pressure balance can establish and 
sonic turbulence remains subdominant, leaving the possi- 
bility for a direct cascade of ion-cyclotron modes. 



As already mentioned, for oblique propagation, the sys- 
tem is driven on the z component of the velocity or the 
magnetic field, thus perpendicularly to the plane defined 
by the ambient field and the direction of propagation. The 
injection rate is the same as the injection on each compo- 
nent in the case of parallel propagation, and thus globally 
reduced by a factor two. Two typical angles are consid- 
ered, 9 = 45" to exemplify a generic oblique regime, and 
9 = 80" that corresponds to a quasi-transverse dynamics. 

4-.1. The case o/45" propagation angle 

With a kinetic driving, the same computational domains 
and forcing scales as in Sections 3.1 and 3.2 were consid- 
ered, while in the case of a magnetic driving, only the 
largest domain was considered. A resolution of 512 mesh 
points was used in all cases, with = pLy = Ky = 2 • 10~^, 
pi^ = K,z — 10~^ and At = 10""* (large-scale forcing), 
fj,^ = ^y — Ky = 5-10^'^, — Hz — 10^^ and At — 5-10^^ 
(smaller-scale forcing) . 

A main observation common to all these simulations 
is the absence of a significant energy transfer to the 
largest scales, but rather the establishment of a statisti- 
cally steady state, with a total energy saturating at about 
the same level in all the runs, after about a thousand 
time units. In the case of the kinetic driving at large 
scales (run E), a pressure balance nevertheless establishes, 
perturbed by localized events associated with the forma- 
tion of cusps on the longitudinal velocity profile (fig. [T5)) . 
These peculiar structures are reminiscent of the interme- 
diate shocks that develop i n the Cohen-Kulsrud equation 
( Cohen and KulsrudL HqtI . The pressure balance dete- 
riorates when the driving is at smaller scales, especially 
when it is of magnetic type. 

At the level of the dissipations, the various types of driv- 
ing do not induce qualitative differences. It is noticeable 
that the dissipation = J {fJ-zidxVz)^ + i^z{dxbz)'^)dx af- 
fecting the velocity and magnetic field components in the 
transverse direction is smaller by a factor of order 10'* 
than the dissipation IJpiano = / {lJ-x{dxVxf + l^y{dxVyY + 
Kyidxby)'^)dx of the components in the (Bo,x) plane. 

A more detailed understanding of the dynamics is pro- 
vided by the frequency analysis displayed in fig. |9l With a 
large-scale kinetic forcing (run E: panels g, h, i), three dis- 
tinct frequencies corresponding to the three linear eigen- 
modes with comparable amplitudes, are identified on the 
k — 1 mode. At smaller scales {k = 4 and k — 25), in con- 
trast, a peak near the intermediate frequency is dominant 
and strongly broadened, with most of the energy contained 
in the component parallel to the direction of propagation, 
consistent with the properties of the corresponding eigen- 
modes displayed in fig. [T] (left). This confirms the com- 
pressible nature of the turbulence that develops at small 
scales. 

When kinetic energy is injected at smaller scales (run 
F: panels q, r), the three eigenfrequencies are still visi- 
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ble at fc = 4, with a much more intense contribution of 
the whistler modes when compared with the large-scale 
driving. The presence of whistler at large or intermediate 
scales which is even more conspicuous in the case of mag- 
netic driving (run G; panels s, t), is expected to be at the 
origin of the disruption of the pressure balance, as also 
observed in the case of parallel propagation. At smaller 
scale (fc ~ 25), the intermediate mode strongly dominates, 
as for the large-scale driving. In other words, in oblique 
directions the presence of dispersion makes the dynamics 
fully compressible, in contrast with parallel propagation 
where an incompressible Alfvenic turbulence can be iso- 
lated when concentrating on the transverse components of 
the fields. 

4.2. The case of 80° propagation angle 

In this case, the simulation was performed with a large- 
scale kinetic driving (run H), in the same conditions as in 
Section 4.1. A resolution of 512 mesh points was used, with 
fJ'x = fJ'y — Ky — 10~^, Hz = = 10^^ and At = 10~^. 
Compared with the 45° case, the energy saturates at a 
much later time {t 40000). This long time interval is 
required because of the presence of an early inverse cas- 
cade that leads to a dominance of the n = 1 mode, which 
nevertheless saturates before the end of the simulation. A 
pressure-balanced state establishes very rapidly (fig. [HI 
left). It is nevertheless perturbed by dispersive shocks, 
which occur very frequently and are especially visible on 
the velocity component , as displayed in fig. [T3] (right) 
at t = 2472. The spectra of the transverse velocity and 
magnetic fields are displayed in fig. 1151 averaged over the 
time intervals t = 2465-2477 (left) and t = 10595-10565 
(right). Left panel represents a typical situation in which a 
dispersive shock is crossing the simulation box, its spectral 
signature being conspicuous as a flat zone at intermediate 
scales. Right panel represents instead the same spectra 
at a moment when no shocks are present. In such a case, 
the velocity spectrum is steeper than the magnetic one, in 
contrast with the case of parallel propagation where the 
opposite takes place. This is consistent with the domina- 
tion of whistler modes at intermediate scales, as discussed 
below. 

Inspection of fig. [S] (panels j, k, 1) provides further 
insights on the dynamics. Compared with the 45° case 
(panels g, h, i), much more energy is contained at low 
frequency in all the spatial Fourier modes. The character- 
istics of the various field at large (fc = 1) and intermediate 
(fc = 4) scales where Vy ~ Vz ^ Vx,by,bz, indicate that 
these modes can be viewed as ion-cyclotron waves (fig. 
[1]), which for this angle of propagation identifies with the 
branch of kinetic Alfven waves. At smaller scales (fc = 25, 
panel 1) Vy ^ f z , probably because of the choice ^y ^ fiz- 
The oscillations in the dispersive shock displayed in fig. [T3] 
(right) can be identified as the whistler modes associated 
with the broadened peak visible on the field for fc = 4 
displayed in panel k. These high-frequency modes are in- 
deed strongly dominated by the Vx component at these 



scales, which is consistent with the properties of the linear 
whistler waves. As a consequence, they are not expected to 
significantly disrupt the pressure balance, in contrast with 
the 45° case with small-scale forcing. In the latter simu- 
lation, whistler modes are indeed present at scales where 
they involve a significant magnetic contribution. At fc = 4, 
a weak peak corresponding to the intermediate frequency 
is also present. The corresponding mode dominates over 
the whistler wave at fc = 25, where it reaches an ampli- 
tude comparable to the low-frequency components (panel 
1). The overall dynamics is then characterized by a com- 
plex distribution of energy among the different waves, in 
contrast with the 45° case where the dynamics is mostly 
governed by the intermediate modes. 

Examination of the dissipations shows no difference be- 
tween the 45° and the 80° cases for what concerns Upiane, 
while Dz is in the mean larger by a factor 2 and dis- 
plays stronger fluctuations for propagation at 80°. The 
presence of ion-cyclotron waves for quasi-transverse prop- 
agation without a major influence on the dissipation 
(which is indeed comparable to that observed at 45° where 
such modes are absent) suggests that they are not able to 
cascade efficiently to the dissipation scales. 

5. Beyond the Hall-MHD description 

Although significantly richer than classical MHD, 
HMHD cannot capture important dynamical properties 
of space plasmas originating from their quasi-collisionless 
character. Although such media should require a fully 
kinetic description and thus enormous computational re- 
sources, it may be of interest to deal with fluid models that 
extend MHD by retaining main kinetic effects such as Lan- 
dau damping and also finite Larmo r radius (FLR) correc- 
tions that, as stressed for example in lServidio et al.l (|2007r ). 
are expected to be comparable to the Hall effect within the 
magnetopause. In spite of their scmi-phenomenological 
nature, these so-called Landau fluid models correctly re- 
produce the kinetic theory in the low-frequency limit. 
First introduced in the context of large-scale MHD 



( Snyder et al.l . ll997f) . these models were then extended by 



retaining Hall effect and FLR corrections in order to de- 
scribe quasi-transverse ionic scales obeying the gyrokinetic 
scaling and thus associated with low frequencies. Never- 
theless, in contrast with the gyrokinetic theory, the fast 
magnetosonic modes arc not averaged out in the Landau 
fluid approach and their large-scale dynamics is accurately 
reprodu ced. A compr ehensive derivation of the model is 



given in lPassot and Sulem (.2007.) , and its reduction to o ne 



space dimension can be found in iBoreoeno et aL (2007). 

Preliminary simulations of the Landau fluid model, re- 
taining Landau damping of the ions and the electrons, 
together with ion finite Larmor radius corrections, were 
performed in a domain of extension IBtt for a propagation 
angle of 80°. The Landau damping being relatively weak in 
quasi-transverse directions, we resorted to add a magnetic 
diffusivity (77 — 0.02) together with a fc^-hyperviscosity 
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Vfi — 10^^ in a simulation with a resolution of 256 collo- 
cation points. The system is driven by a random forcing 
acting on the transverse components of the velocity field, 
identical to that used in the HMHD simulation presented 
in section 4.2. Nevertheless, in order to control the growth 
of the total energy that is dominated by the global heat- 
ing of the plasma, the forcing is monitored in order to 
maintain the kinetic energy between prescribed upper and 
lower bounds. 

As in the HMHD corresponding simulation, we ob- 
serve the development of pressure-balanced structures 
with the formation of small-scale fluctuations associated 
with whistler modes. Compressible modes are neverthe- 
less more damped than in HMHD simulations due to Lan- 
dau damping. An interesting issue that the Landau fiuid 
can address concerns the evolution of the plasma temper- 
atures. Figure [TBI displays the typical evolution of the ion 
temperatures of a plasma with ion parallel beta /3||j equal 
to 0.6 (left) and 1.2 (right). We observe in both simula- 
tions that the anisotropy progressively grows up, leading to 
dominant transverse ion temperatures, making the plasma 
potentially unstable to ion mirror instabilities. For the 
electrons, the parallel temperature (not shown) increases, 
possibly leading to an electron firehose instability. Later 
on, we observe abrupt variations of parallel ion temper- 
atures, resulting in a reduction of the anisotropy. These 
events are associated with the formation of quasi-singular 
structures on the velocity components, that can be shown 
to affect more the parallel than the perpendicular tem- 
perature. The above temperature variations was checked 
to be due to mirror instabilities, that is saturated thanks 
to the small-scale FLR corrections retained by the model. 
Differently, a phenomenological relaxation of temperature 
instabilitie s, in the form of effe ctive collisions, was imple- 
mented in [Sh arma et alj I (l2006l) . Such a model should be 
useful in situations where ion-cyclotron or oblique firehose 
instabilities are excited, as these instabilities cannot be 
described with fluid equations. 

The Landau fluid simulations thus reproduce constrain- 
ing effects due to temperature anisotropy instabilities on 
the plasma parameter s, also observed on solar wind data 



(jHellinger et al.l . l2006f ) . As the plasma anisotropy relaxes 



under the effects of temperature instabilities, the system 
enters a phase of slow dynamics. This regime is however 
not expected to persist, as new temperature anisotropy 
can develop, making the above scenario to repeat. 

6. Conclusion 

Although limited to one space dimension, the present 
study reveals specific aspects of the turbulent dynamics of 
magnetized plasmas at scales comparable to the ion iner- 
tial length, in a regime where the transverse components 
of the velocity or the magnetic fields are randomly driven. 
Special attention was paid to the distribution of the en- 
ergy among the different MHD waves that can be clearly 
identified from their linear dispersion relation, in spite of a 



possible small shift in the temporal spectrum of the com- 
puted fields, in situations where the presence of large-scale 
coherent structures can be viewed as performing a renor- 
malization of the ambient parameters. 

A main observation in the case of parallel propagation is 
the contrast between a large-scale forcing for which the en- 
ergy is almost entirely transferred to magnetosonic modes 
with nevertheless a non- negligible transfer to larger scales, 
and the regime where the driving takes place at scales 
where dispersion is more efficient, for which, provided the 
driving is of kinetic type, a direct Alfvenic transfer es- 
tablishes, a result qualitatively co nsistent with the weak- 
turbulence analysis performed by lYoon and Fang 1200811 
in the context of the Vlasov equation. In oblique direc- 
tions, the combined role of dispersion and compressibil- 
ity leads to a turbulence dominated by the intermediate 
modes, with an increasing contribution of low frequency 
kinetic Alfven waves, together with a faster establishment 
of total pressure balance, as the propagation angle is in- 
creased. Furthermore, the Landau fiuid model shades a 
light on the development of pressure anisotropy resulting 
in recurrent instabilities. 

To conclude, we would like to stress the complexity 
of the turbulence problem in magnetized plasmas, even 
within the strongly simplified description provided by one- 
dimensional HMHD. The usual picture of incrtial ranges 
where energy "cascades" progressively from scale to scale 
at a constant rate turns out to be strongly affected by the 
predominance of structures and by a strongly fluctuating 
transfer, making possibly questionable the usual concepts 
of the classical turbulence theory, and leading to a dynam- 
ics that turns out to be significantly less universal. 
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Figure 1: Wavenumber dependency of the linear HMHD eigenfre- 
quencies (top) and eigenmode components W'^^ , W'-^^ , W'^^ (from 
middle to bottom), when prescribing p = 1, for propagation angles 
e = 45° (left) and 6 = 80° (right): magnetic field by (thick solid line) 
and bz (thick dashed line); velocity components Vy (thin solid line) 
and Vz (thin dashed line); velocity component Vx (dotted line). 




Figure 3: Left; profile of the the magnetic pressure |b^p/2 (solid 
line) and of the thermal pressure ^fP (dashed line) , together with 
their half-sum (dotted line) showing pressure balance at time t = 
20000 for the run A. Right: transverse velocity field intensity |vxP 
at the same time. 




Figure 4: Transverse velocity field intensity (IvxP) (left) and indi- 
vidual components Vy , Vz at time t = 18200 for run A, after filtering 
the mode of index n = 1. 
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Figure 2: Left: time evolution of the total energy and of the trans- 
verse kinetic and magnetic contributions for run A with kinetic forc- 
ing discussed in Section 3.1. Right; time evolution of the magnetic 
energy |b^p -|- \hz\'^ for the Fourier modes of index 1 (solid line), 
2(dashed), 3(dotted), 4 (dash-dotted), for the same simulation. 



Figure 5: Kinetic |vyp + |v2p (thick line) and magnetic Ib^p-l-jb^p 
(thin line) spectra averaged over the time interval t = 18200 — 18250 
(left), and t = 19500 - 19550 (right) for run A. 
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Figure 6: Left: time evolution of the parallel dissipation Dy for 
the kinetic (solid line) and magnetic (dashed line) drivings (runs A 
and B respectively) described in Section 3.1. Right: same for the 
perpendicular dissipation D±. 
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Figure 7: Time evolution of the parallel dissipation (dark line) and 
of the rate of energy transfer S from the Alfven to magnetosonic 
waves (light line) for run A. Left: instantaneous quantities. Right: 
time-averaged quantities over Atave = 1000. 
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Figure 12: Left: time evolution of the parallel dissipation D|| for 
the kinetic (run C, solid line) and magnetic (run D, dashed line) 
drivings in the conditions described in Section 3.2. Right: same for 
the perpendicular dissipation Dj_ (in lin-log scale). 
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Figure 8: Left: kinetic Ivyp + lv^p (thick line) and magnetic |bj;p + 
|bzp (thin line) spectra averaged over the time interval t = 28000 — 
28050 for the kinetically driven run of Section 3.1 with monitored 
forcing. Right: transverse velocity field intensity (|vxP) at time 
t = 28000, in the same conditions. 




Figure 13: Left: profile of the the magnetic |bxP/2 (solid line) and 
thermal 13 p'^ /-y (dashed line) pressures, together with their half-sum 
(dotted line) at time t = 422, for the run E with large-scale kinetic 
driving described in Section 4.1. Right: velocity component Vx at 
the same time. 
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Figure 10: Time evolution of the total kinetic and magnetic ener- 
gies for kinetic (left) or magnetic (right) drivings (runs C and D 
respectively) in the conditions described in Section 3.2. 
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Figure 14: Left: profile of the the magnetic |bxP/2 (solid line) and 
thermal ,Sp'^/7 (dashed line) pressures, together with their half-sum 
(dotted line) at time t = 2472 for run H with large-scale kinetic 
driving described in Section 4.2. Right: velocity component Vx at 
the same time. 




Figure 11: Profile of |vxp at t = 4000 for kinetic (run C, left) and 
magnetic (run D, right) driving for the run described in Section 3.2. 




Figure 15: Spectra |vzp (thick line) and Ib^p (thin line) of the 
transverse fields, averaged over the time intervals t = 2465—24772477 
(left), and t = 10595 - 10565 (right) for run H. 
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Figure 16: Time evolution of the parallel (grey line) and transverse 
(black line) temperatures for simulations with /3||i = 0.6 (left) and 
= 1.2 (right) for the runs described in Section 5. 
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